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We summarize recent progress in soft QCD modeling based on the set of Dyson-Schwinger equa- 
tions truncated to ladder-rainbow level. This covariant approach to hadron physics accommodates 
quark confinement and implements the QCD one-loop renormalization group behavior. We compare 
the dressed quark propagator, pseudoscalar and vector meson masses as a function of quark mass, 
and the p — > tttv coupling to recent lattice-QCD data. The error in the Gell-Mann-Oakes-Renner 
relation with increasing quark mass is quantified by comparison to the exact pseudoscalar mass 
relation as evaluated within the ladder-rainbow Dyson-Schwinger model. 
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I. INTRODUCTION 

The study of light-quark pseudoscalar and vector 
mesons is an important tool for understanding how QCD 
works in the non-perturbative regime. The pseudoscalars 
are important because they are the lightest observed 
hadrons and are the Goldstone bosons associated with 
dynamical chiral symmetry breaking. The ground state 
vector mesons are important because, as the lowest spin 
excitations of the pseudoscalars, they relate closely to 
hadronic qq modes that are electromagnetically excited. 

We use a Poincare covariant model defined within the 
framework of the Dyson-Schwinger equations [DSEs] of 
QCD; these form an excellent tool to study nonperturba- 
tive aspects of hadron properties |]j . It is straightforward 
to implement the correct one-loop renormalization group 
behavior of QCD Q , and obtain agreement with pertur- 
bation theory in the perturbative region. Provided that 
the relevant Ward-Takahashi identities are preserved in 
the truncation of the DSEs, the corresponding currents 
are conserved. Axial current conservation induces the 
Goldstone nature of the pions and kaons ; electromag- 
netic current conservation produces the correct electric 
charge of the mesons without fine-tuning. These proper- 
ties are implemented here within the rainbow truncation 
of the DSE for the dressed quark propagators together 
with the ladder approximation for the Bethe-Salpeter 
equation [BSE] for meson bound states. 

The model we use has two infrared parameters 
which specify the momentum distribution and strength 
of the ladder-rainbow kernel at a low scale necessary to 
generate an empirically acceptable amount of dynamical 
chiral symmetry breaking [3, |6( as measured by the chi- 
ral condensate. As a corollary, the strong dressing of the 
quark propagator shifts the mass pole significantly away 
from the real timelike p 2 axis. The produced bound state 
mesons do not have a qq decay width and, in this sense, 
the present model implements quark confinement. The 



absence of a real mass pole for dressed quark and gluon 
propagators has been studied and found to be a sufficient, 
but not necessary, condition for confinement [ij, |3, la, |9( . 
The model provides an efficient description of the masses 
and decay constants of the light-quark pseudoscalar and 
vector mesons 0, Q-the elastic charge form factors 
F 7T (Q 2 ) and Fk(Q 2 ) [l3 and the electroweak transition 
form factors of the pseudoscalars and vectors [Til Il2| . 



II. DYSON-SCHWINGER EQUATIONS 



The dressed quark propagator S(p) is the solution to 
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+Zi J g'D^(p-q)-^S(q)Tl(q,p) (1) 

where Dy, v (k) is the renormalized dressed-gluon prop- 
agator, Tl(q,p) is the renormalized dressed quark-gluon 
vertex. The notation = J A d 4 g/(27r) 4 denotes a trans- 
lationally invariant regularization of the integral with 
mass-scale A. The solution is renormalized according to 
Sip) -1 — vy-p-{-m(fJb) at a sufficiently large spacelike /i 2 , 
with m(/x) the renormalized quark mass at the scale \i. 
The renormalization constants Z2 and Z^ depend on the 
renormalization mass-scale /j, and on the regularization 
mass-scale A. The limit A — > 00 is to be taken at the end 
of all calculations. 

The BSE for a ab meson is 



T ab (p + ,p_) = / K(p,q;P)S a (q + )T ab ( q+ , q _)S b (q-), 

Jq 

where K is the renormalized qq scattering kernel that is 
irreducible with respect to a pair of qq lines. The quark 
momenta are q± ; the meson momentum is P = q+ — q- 
and satisfies P 2 = —m 2 . The relative momentum q is 
introduced by q + = q + r\P and </_ = q— (1 — rj)P where 
•q is the momentum partitioning parameter. Physical ob- 
servables should not depend on 77 and this provides a 
convenient check on numerical methods. We employ the 
model that has been developed recently for an efficient 
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FIG. 1: DSE solution plll3f for qua rk p ropagator amplitudes 
compared to recent lattice data (l4lll5j|. 



description of the masses and decay constants of the light 
pseudoscalar and vector mesons 0, @| • This consists of 
the rainbow truncation of the DSE for the quark propaga- 
tor and the ladder truncation of the BSE for the pion and 
kaon amplitudes. The required effective qq interaction is 
constrained by perturbative QCD in the ultraviolet and 
has a phenomenological infrared behavior. In particular, 
the rainbow truncation of the quark DSE, Eq. Q), and 
the ladder truncation of the BSE, Eq. @, are 



Z ig 2 D^(k)Ti(q,p) 



4™ off (fc 2 )I?^ (fc)7^ 



and 



K(p,q;P) 



-Ana cS {k 2 )D^ c (k) 



(3) 



(4) 



where D l *° c (k = p — q) is the free gluon propagator in 
Landau gauge. These two truncations are consistent in 
the sense that the combination produces vector and axial- 
vector vertices satisfying the respective WTIs. In the ax- 
ial case, this ensures that in the chiral limit the ground 
state pseudoscalar mesons are the massless Goldstone 
bosons associated with chiral symmetry breaking 0, • 
In the vector case, this ensures electromagnetic current 
conservation. The "effective coupling" a c ff(fc 2 ) defines 
the model. The ultraviolet behavior is chosen to be that 
of the QCD running coupling a(k 2 ); the ladder-rainbow 
truncation then generates the correct perturbative QCD 
structure of the DSE-BSE system of equations. The phe- 
nomenological infrared form of a e g (k 2 ) is chosen so that 
the DSE kernel contains sufficient infrared enhancement 
to produce an empirically acceptable amount of dynami- 
cal chiral symmetry breaking as represented by the chiral 
condensate [iH ]. 

We employ the Ansatz found to be successful in earlier 



TABLE I: The pseudoscalar observables that define the 
present ladder-rainbow DSE-BSE model, adapted from 
Refs. 0,0. 
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first term implements the strong infrared enhancement 
in the region < k 2 < 1 GeV 2 required for sufficient dy- 
namical chiral symmetry breaking. The second term 
serves to preserve the one-loop renormalization group 
behavior of QCD. We use m t = 0.5GeV, r = c 2 — 1, 
Nf = 4, and we take Aq C d = 0.234 GeV. The renormal- 
ization scale is chosen to be p — 19 GeV which is well into 
the domain where one-loop perturbative behavior is ap- 
propriate 0,3- The remaining parameters, u = 0.4 GeV 
and D = 0.93 GeV 2 along with the quark masses, are fit- 
ted to give a good description of (qq), rn^/ K and / w as 
shown in Tabled The subsequent values for fx and the 
masses and decay constants of the vector mesons p, <f>, K* 
are found to be within 10% of the experimental data p| . 

In Fig. n we compare the DSE model propagator am- 
plitudes defined by S(p) = Z(p 2 )[i^)+ Mip 2 )}- 1 with 
the most recent results in lattice QCD using staggered 
fermions in Landau gauge |14l Il5| . These simulations 
were done with the Asqtad improved staggered quark ac- 
tion, which has lattice errors of order C(a 4 ) and 0(a 2 g 2 ). 
Fig. n shows both M(p) and Z(p) obtained with a bare 
lattice mass of ma = 0.036 in lattice units, which corre- 
sponds to a bare mass of 57 MeV in physical units. The 
DSE calculations use a current mass value of 75 MeV 
at /i = 1 GeV to match the lattice mass function around 
3 GeV; this current mass is about 0.6 m s . There is agree- 
ment in the qualitative infrared structure of M(p) and 
Z(p). Since the lattice simulation produces the regulated 
but un-renormalized propagator, the scale of Z(p) is arbi- 
trary and we have rescaled the lattice Z(p) to match the 
DSE solution at 3 GeV. For Z(p) the ladder-rainbow DSE 
model saturates much slower than does the lattice; this 
may signal a deficiency of the bare gluon-quark vertex. A 
recent study of the coupled ghost-gluon-quark DSEs has 
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FIG. 2: DSE calculation of equal flavor vector meson mass 
variation with pseudoscalar meson mass as m q is varied com- 
pared to lattice data from CP-PACS EH and UK-QCD [l7|. 



found that the quark-gluon vertex dressing can produce 
a change of this character in the infrared structure of the 
quark amplitudes 0, |2(| • 

In Fig. |2 we compare the rainbow-ladder DSE model 
with unquenched lattice data for the variation of vector 
and pseudoscalar meson masses with quark current mass 
in the case of equal flavor quarks. The DSE calculation, 
shown by the discrete circles, is limited to the mass range 
where it is reliable. The solid line is a fit to those results 
plus the experimental J/tp — r\ c point. The curvature at 
low mass is consistent with My cx m q and Mp cx ^Jrn q . 
This comparison is consistent with the known properties 
of the DSE results: My is 5% too low for the p and 5% 
too high for the (j>, while the pseudoscalar masses in the 
u-quark and s-quark regions are fit to experiment. (If 
the evident fit is continued to the T vector mass, the 
predicted r\h mass would be 10.0 GeV.) 



III. PSEUDOSCALAR MESON MASS 
RELATION 

As the current quark mass is raised from zero, the ex- 
plicit breaking of chiral symmetry adds mass to the Gold- 
stone boson modes. The way in which the pseudoscalar 
meson mass grows with quark mass is described, at low 
mass, by the GMOR relation. This is 



Mp(m 1 (/j),m 2 (/i)) = [mi(^)+m 2 (/i)] 



(ft 



-0{m 2 ) 



(6) 

in the general case where the two quark flavors are dif- 
ferent. Here (qq)^ = —Z<±N c tr J A So(q) is the chiral con- 
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FIG. 3: Mv(rriQ) and Mp(rriQ) for unequal flavor u — Q 
mesons. The DSE model-exact mass relation reveals the size 
of the correction to the GMOR relation. 



densate at scale /z, the current masses are determined at 
the same scale, and f P is the chiral limit electroweak de- 
cay constant (in the f n = 92.4 MeV convention). For u/d 
quarks the GMOR relation is satisfied to high accuracy 
(within 0.2% in the present DSE model). For current 
masses of the order of m s ~ 120 MeV and above, the 
question of the size of the error in the GMOR relation 
is not settled. An exact mass relation for pseudoscalar 
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FIG. 4: Current quark mass dependence of fp, the elec- 
troweak decay constant for u — Q mesons. [f„ — 131 MeV 
convention.] 

mesons in QCD, applicable for all values of the quark 
masses, has been established p. With allowance for dif- 
ferent quark flavors, it takes the form 

Mp /p(mi,m 2 ) = (mi + m 2 ) Rp (mi,m 2 ) , (7) 
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where Rp is the projection of the meson wave function 
onto 75 at the origin of qq separation and is given by |3| 

,A 

Rp = -iZ 4 N c tr[ 75 S fl (q + )r P (q;P)S h (q^], (8) 

Jq 

with all renormalized quantities taken at the same scale 
fi. The chiral limit of this quantity can be shown to 
be Rp{m x = m 2 = 0) = -(qq)°/f P , and thus the GMOR 
relation follows as a collorary of the exact relation, 
Eq. J2J, at low mass. The origin of this exact mass rela- 
tion is the axial vector Ward-Takahashi relation 

- iP? T 5fl (q; P) = S£{q+) 75 + 75 Sj*(q-) 

-(mi + m 2 )r 5 (g;P) . (9) 

The dressed vertices T^^q; P) and T^(q; P) satisfy inho- 
mogeneous integral equations that have the same kernel, 
the irreducible qq scattering amplitude; the inhomoge- 
neous terms are Z 2 7m and Z 4 75 respectively. Thus both 
vertices have poles corresponding to the pseudoscalar 
meson bound states. [The axial vector poles in T^,^ 
have transverse residues and do not contribute.] The 
exact mass relation, Eq. (JJJ), arises from the equality 
of the pseudoscalar pole residues from both sides of 
Eq. ©. The residue of T 5 (q;P) is -i R P T P (q;P), and 
the residue of T 5ll (q;P) is P M fpTp(q; P), where Tp is 
the pseudoscalar bound state BS amplitude. The ex- 
pression for fpP^ is the same as for Rp in Eq. © ex- 
cept that 75 is replaced by 75 7 M and —iZ 4 is replaced 
by Zij\rl (in the convention where the physical f n is 
92.4 MeV). In the chiral limit, the last term of Eq. © 
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FIG. 5: A measure of the current quark mass dependence of 
the in-hadron condensate for qQ pseudoscalars. 

is not present and the right hand side is not singular. A 
systematic expansion in powers of P^ reveals |3| that in 
leading order the pole in the left hand side must move to 



P 2 = to be cancelled and thus Mp = 0. The analysis 
also provides relationships between the BS amplitude and 
the quark propagator, e.g., f P E p (q;0) = E>o(q 2 ), where 
Tp(g; P) — ij^Ep + • • •. Since the quark mass function 
is proportional to Bo(q 2 ), this latter relation means that 
dynamical chiral symmetry breaking is necessarily ac- 
companied by a massless pseudoscalar bound state. This 
is the familar Goldstone's theorem; notice that the com- 
posite, distributed nature of the pion amplitude requires 
a running quark mass function. 

At small m, both fp and Rp are constant leading to 
the GMOR behavior Mp ~ s/m. The error in this has 
to increase with mass since the heavy quark limiting be- 
havior is I2H fp ~ 1/s/fri and Rp ~ \fm which leads to 
the linear behavior Mp ~ to. 
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FIG. 6: Our DSE results for the symmetric 7*7r — > 7* form 
factor, compared to the pQCD asymptotic 1/Q 2 behavior. 
[Here f n — 131 MeV.] The naive VMD model suggests a 
dipole behavior which is correct only in the infrared. 

The evolution of both vector and pseudoscalar meson 
masses with increasing current quark mass has been stud- 
ied within the DSE model in both the equal and unequal 
flavor cases. An example is shown in Fig. [S] for u — Q 
mesons as a function of vtlq (/1 — 1 GeV) up to the limit 
of accuracy of the calculations. The fit shown is adapted 
from Ref. (22J by evolving the m<g(/i = 19 GeV) used 
there to \i = 1 GeV for ease of comparison with conven- 
tionally quoted values. Thus here we have Mp = a'+ 
/3'y/mQ+ 'j' m Q'i with both masses in GeV, the parame- 
ters are (a', /?', 7') = (0.083, 0.842, 0.880). With the DSE 
model-exact values of f P and (qq)^, the GMOR relation 
is explicitly Mp = 0.00955+ 1.724 m Q (fj. = 1 GeV) and 
this is compared to the exact mass relation in Fig. El For 
the K meson the GMOR error is 4%, at m Q = 0.4 GeV 
the error is 14%, while at the D meson the error is 30%. 
In the D meson region, only ~ 50% of the mass comes 
from the linear term; the heavy quark domain is at higher 
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TABLE II: The coupling constants g v ->pp calculated in the 
DSE model compared to results from a fit to the timelike form 
factor pole plE^ and lattice-QCD [2(J. 



Qv^pp 


Expt 


this work 


pole fit 


lattice-QCD 


Qp — >7T7T 


6.02 


5.14 


5.2 


6-08±; : gg 




4.64 


4.25 


4.3 




9k* + ^,k0tt+ 


4.60 


4.81 


4.1 





mass. 

The exact mass relation Mp(ijiq), Eq. Q, differs from 
the GMOR behavior due to the mass dependence of fp 
and Rp. Instead of the latter quantity, one can de- 
fine (qq) ^ = — fp Rp as an effective "in-hadron" conden- 
sate |21| which allows the exact mass relation to take the 
GMOR-like form Alp fp = (mi + rri 2 ) f , \{qq)"\. From 
this relation we extract the quark mass dependence of 
(m)" at fi = 1 GeV. The low mass DSE results for M P 
and fp are fitted to forms that respect the heavy quark 
limits, are consistent with D and B meson masses, and 
are consistent with lattice results j2^| for fp of the D 
and B mesons. We adapt the fp(Mp) fit from Ref. [2l| 
to produce fp(rriQ) and to accommodate the DSE re- 
sults at low mass. The result is displayed in Fig. 0] 
The fit (in the f n = 131 MeV convention) is fp = N/D, 
where N = a + bniQ and D = 1.0 + ctuq + dm,Q, with 
(a,b,c,d)= (0.017,0.068,0.649,0.391) when m Q is in 
GeV. The resulting estimate of the mass dependence for 
(qq)Jf at \i = 1 GeV is shown in Fig. [S] and indicates 
~ 15% increase over the chiral limit value for an tiiq 
relevant to K, while for the D meson the in-hadron con- 
densate is about 70% enhanced. 



IV. VECTOR MESON STRONG DECAYS 

Since this DSE model describes the elastic charge form 
factors of the pseudoscalars very well 0] m impulse 
approximation, the strong decays of the vector mesons 
should be well-described without parameter adjustment. 
In impulse approximation, the amplitude for the decay of 
p with 4-momentum Q = pi + P2 to 7T7T with 4-momenta 
Pi,P2 is given by [l| 

,A 

2PJ g P ^ = V2N c tr J S(q) Tp(q 7 g+) S(q+) 

xr M ( g+ ,g_)%_)f P (g_,g), (10) 

where no distinction is made between the u/d quarks, 
P = (pi - p 2 )/2, q = k + P/2 and q ± = k - P/2± Q/2. 
The appropriate generalizations of the flavor structure 
appropriate to <f> — > KK and K* — > Kir are straightfor- 
ward |l3||. In Eq. (|TU|l . the component of P M transverse 
to Qfj, is indicated on the left hand side to cover the case 
of unequal decay products. 

The results shown in Table ITD are within 5-10% of ex- 
periment with the error being larger if the vector meson 



is lighter. As an independent check, coupling constants 
are also extracted from the timelike electroweak form fac- 
tors near the vector meson poles [3, . The agreement 
is encouraging considering that with eight independent 
covariants for the vector BS amplitude and four each for 
the pseudoscalars, there are 128 distinct quark loop in- 
tegrals for each physical decay. Also shown in Table |n] 
is a recent lattice-QCD result [26J for p — > 7T7t from the 
UKQCD Collaboration. Although the lattice data is at 
m^/m p — 0.578, which corresponds to the s-quark mass, 
and thus no physical decay of the p can take place, the 
amplitude (p|7T7r) is accessible through study of state mix- 
ing on the lattice. 

Since the width of the p is almost 20% of its mass 
while the widths of the <f> and K* are significantly less 
important, we expect the ladder approximation for the 
BSE kernel (which omits the strong channels 7T7t, KK 
and Kir respectively) to be less accurate for the p than 
for the (j> and K*. Accordingly we speculate that this is 
largely the reason why the result for gp^-n-K in Table |n] 
deviates from experiment twice as much (15%) as do the 
other decay constants. 



V. PQCD LIMIT OF FORM FACTORS 

Besides the soft physical characteristics of light 
mesons, the present DSE model should also reproduce 
perturbative QCD limits. This has been checked for the 
uv behavior of the quark mass function M(p 2 ); both the 
leading log behavior away from the chiral limit, and the 
coefficient of the leading 1 /p 2 behavior in the chiral limit 
reproduce the exact 1-loop results of QCD 0- A more 
difficult task is to test the asymptotic behavior of meson 
form factors against pQCD predictions. This is compli- 
cated by the fact that covariant ladder-rainbow calcula- 
tions that link the dressed quark propagator, the BSE, 
and the impulse approximation for form factors have only 
been carried out in Euclidean metric for practical reasons. 
The mass-shell constraint for mesons then requires an 
analytic continuation which entails complex quark mo- 
menta in loop integrals. This greatly hinders the asymp- 
totic analysis. 

A case that is free of these difficulties is the symmet- 
ric 7*7r — ► 7* transition where the photons are taken 
to have equal virtuality Q 2 and there is only one mass- 
shell constraint. Since m 2 . is negligible compared to all 
other scales in the problem, all involved quark momenta 
are essentially real and spacelike. In Fig. we show the 
result 0| of the present DSE model compared to the 
pQCD asymptotic behavior (2?J obtained from the light- 
cone operator product expansion. (In this case, log cor- 
rections occur at sub- leading order.) The numerically 
generated asymptotic behavior of the DSE-based model 
reproduces the pQCD limit as it must. By about 2 GeV 2 
the dressing of the photon vertices becomes negligible; 
however the 3-point function does not become an effective 
2-point function (thereby generating the required power 
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of fa) until about 15-20 GeV 2 HJ. Such a high scale for 
the onset of pQCD behavior is consistent with an earlier 
observation |23 hi a DSE-based model study of F^(Q 2 ). 

VI. DISCUSSION 

Recent reviews 0, put this model in a wider per- 
spective and compile results for both meson and baryon 
physics, an analysis how quark confinement is manifest in 
solutions of the DSEs, and both finite temperature and 
finite density extensions. The question of the relevance 
and accuracy of the ladder-rainbow truncation has also 
received some attention; it has been shown to be partic- 
ularly suitable for the flavor octet pseudoscalar mesons 
since the next-order contributions to the BSE kernel, in a 
quark-gluon skeleton graph expansion, have a significant 



amount of cancellation between repulsive and attractive 
corrections [jlH EI ■ The preservation of the axial vec- 
tor WTI is what makes the pseudoscalar meson sector a 
robust and ideal base for parameter fixing; the rainbow- 
ladder truncation may be used as a convenience in that 
sector. It is hoped that future interplay between lattice 
simulations and continuum modeling will increase our un- 
derstanding of QCD for hadron physics. 
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